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Abstract 



Suppose one buys two very similar stocks and is curious about how much, after some time T, one 
of them will contribute to the overall asset, expecting, of course, that it should be around 1/2 of the 
sum. Here we examine this question within the classical Black and Scholes (BS) model, focusing on 
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the evolution of the probability density function P{w) of a random variable w = Uj, /{aj, + Qj. ) 
where a^ and a^ are the values of two (either European- or the Asian-style) options produced by 
two absolutely identical BS stochastic equations. We show that within the realm of the BS model the 
behavior of P{w) is surprisingly different from common-sense-based expectations. For the European- 
Q^l style options P{'w) always undergoes a transition, (when T approaches a certain threshold value), from 

a unimodal to a bimodal form with the most probable values being close to and 1, and, strikingly, w = 

p^ I 1/2 being the least probable value. This signifies that the symmetry between two options spontaneously 

^ ■ breaks and just one of them completely dominates the sum. For path-dependent Asian-style options we 

^^ I observe the same anomalous behavior, but only for a certain range of parameters. Outside of this range, 

>0 ■ P(w) is always a bell-shaped function with a maximum at w = 1/2. 

O ■ Introduction 

O' 



In finance, the style of an option is a general term denoting the class into which the option falls, usually 
defined by the dates on which the option may be exercised. An Asian option is an option where the payoff is 
not determined by the underlying price at maturity, contrary to European or American-style options, but by 
C^ I the average underlying price over some pre-set period of time. Such average-value options are commonly 

traded on some Asian (but also Western) markets being somehow advantageous over European ones since 
the risk of asset price manipulation near the maturity date is reduced due to their path dependence. 

In classical Black-Scholes settings (Black and Scholes 1973, Merton 1973), the underlying asset on which 
the Asian option is based is equal to a stock with price St which follows the so-called geometric Brownian 
motion. In other words, it is defined as the strong solution of the following hnear stochastic differential 
equation: 

dSt = uj Stdt + a StdBt, (1) 

where a is the volatility of St, Bt is a standard one-dimensional Brownian motion and a; is a constant 
dependent on the nature of the asset: it could be a stock, a currency, a commodity and etc. For example, if 
St is a stock paying a dividend at the continuous rate 5, one has u) = r — 5, where r is the risk-free rate of 
interest. 



Resorting to Ito calculus, one solves ^ to find 



^2 



St = So exp ( — — /xt + aBt j , (2) 

where Sq = St=(j is the initial price and /i = 1 — 2uj/a'^ is a constant, which may be positive, equal to 
zero or negative. Further on, for continuous averaging with equal time-independent weights, the random 
variable of interest - the value of an Asian option - is defined as the following functional of a Brownian 
trajectory (Yor 1992, Matsumoto and Yor 2005): 



S f'^ f (J^ \ 

At= J, _°^ r, T= dtexpi — jfit + aBtj 



(3) 



where T is the maturity date and Tq is the time moment when one starts to monitor the evolution of St- 
Without lack of generality, we set Tq = in what follows. 

In this paper we pose the following question which is, we believe, of a considerable conceptual interest: 
Suppose that one has not a single equation ^ but, say, two such stochastic equations, having the same 
volatility, the same u, the same initial price 5*0 and random noise terms, dBl and dB^ , which have 
identic distributions. These two equations produce two identical European-style assets, S\ and S^ , which, 
in turn, generate two identical Asian-style random variables, Aj,' and Aj^. What can be said about the 
distribution functions P{W) and P{W) of random variables: 

c(i) 
W = (4) 
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and 
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These random variables describe a realization-dependent relative weight of one asset in the sum of two 
assets. Noticing that S^ {Aj, ) and S\ {Aj, ) are obviously equal to each other on average, and moreover, 
that all their higher moments are equal, one might be tempted to say that the distribution functions P{W) 
and P{yV) should be bell-shaped functions with a maximum at 1/2. They may broaden (initially both are 
delta-functions) with growth of the maturity T, but still W = 1/2 and W = 1/2 should remain the most 
probable values. 

Curiously, within the realm of the Black-Scholes model most often this is not the case. We set out to show 
here that the behavior of P{W) and PiyV) is, in general, surprisingly different from these common-sense- 
based expectations. 



1 Competition of two uncorrelated European-style variables. 

Note first that W (and hence, P{W)) is independent of jjl, since the factors exjp{—a^iJ.T/2) in the nominator 
and the denominator cancel each other. 

For uncorrelated increments dB^ and dB^ , the distribution function P{W) can be calculated exactly (see 
Appendix lAl): 
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Figure 1: Probability density P(W) in ^ for different values of the effective maturity a^: Dotted curve 
corresponds to ay = 0.25, dashed - to ay = 0.5, while dot-dashed one - to qt = 0.7. 



where ut = a^T jl is an effective maturity. This is the so-called logit-normal distribution, i.e., a probability 
distribution of a random variable whose logit has a normal distribution (Johnson et al 1994). 

Despite a relatively simple form, P{W) in ^ contains a surprise: it shows a markedly different behavior for 
ar less or greater than a^ = 1/2 (see Fig.[T|l. For a^ < 1/2 the distribution has a maximum dXW = 1/2, 
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which means that at early stages both Sj, and Sj, contribute proportionally; it is thus most likely that 
each variable defines just one half of the sum. However, when ar exceeds 1/2, P{W^ changes its shape 
from a unimodal to a bimodal, M-shaped form with P{W = 0) = P{W = 1) = and maximal values 
progressively closer to and 1 as ay — >■ oo. Strikingly, W = 1/2 is now the least probable value. It means 
that at sufficiently large maturities the symmetry breaks and one of the variables completely dominates the 
sum, while the second one becomes a complete loser. Certainly, this is not the behavior one may expect on 
intuitive grounds. 



2 Competition of two uncorrelated Asian-style variables. 



Let ^(r) denote the distribution function of r variables in ^. Then, for two uncorrelated Asian-style 
variables the probability density P(>V) can be written as (see Appendix JB]): 



P(W) 
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Clearly, P{y^^ is symmetric around W = 1/2. The question is whether W = 1/2 is always the maximum 
of the distribution? 

The form of the distribution function ^(r) was discussed in the Uterature on mathematical finance [see, e.g., 
(Yor 1992, Geman and Yor 1993, Dufresne 2004)]. In addition, such r variables appear in different domains 
of probability theory: they are the continuous counterparts of the so-called Kesten variables (Kesten 1973) 
that play an important role in multiplicative stochastic processes and in the renewal theory for products of 
random matrices. In the physical literature, variables r emerge in different contexts related to transport in 
disordered media. In particular, r defines a resistance of a finite interval of length T offered to a passage 
(Redner 2001) of particles diffusing in presence of a random, time-independent Gaussian force with average 



value o"^/z/2 and variance a"^ (Kesten et al 1975, Sinai 1982, Derrida and Pomeau 1982). Consequently, 
inverse moments of r define the moments of stationary currents through the interval boundaries. Within 
this context, moments of r (Burlatsky et al 1992, Oshanin et al 1993a, Oshanin et al 1993b, Monthus and 
Comtet 1994, Comtet et al 1998) and the distribution function *(r) (Monthus and Comtet 1994, Comtet et 
al 1998) were also calculated exactly. 



For arbitrary /x, ^{t) is determined by (Monthus and Comtet 1994, Comtet et al 1998): 

^(r) = ^con(r) + ^di.(r), 
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with 
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where r' = (T^t/2, LZ{x) are generalized Laguerre polynomials, r(x) is the Gamma function, and W^ 



p,'^\ 



are Whittaker functions (Abramowitz and Stegun 1972). For /u > the distribution function in ^ consists 
of discrete ["^dis{T)] and continuous branches ["^con{T)]^ while for ^ < the distribution is determined by 
the continuous branch only. Hence, unlike the distribution of W, P{W) will depend on the sign (and value) 
of fi. Consequently, we will consider the cases of positive and negative fi separately. 
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0. 



In this marginal case r defines an inverse probability current in a finite Sinai chain (Burlatsky et al 1992, 
Oshanin et al 1993a, Oshanin et al 1993b). Occasionally, for this case the distribution P(W) has been 
already calculated in (Oshanin and Redner 2009), which studied the probability that a partially melted 
heteropolymer at the melting temperature will either denaturate completely or return back to a native helix 
state. Here W defines the so-called splitting probability (Redner 2001) - the probability that the boundary 
between the helix and coil phases (which performs Sinai-type diffusion (Sinai 1982)) will first hit one of the 
extremities of the chain without having ever reached the second extremity. 



Adapting to our notations the result of (Oshanin and Redner 2009), we have 
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, rj = arcsinh 
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In Fig. [2] we plot P{yV) for several values of an effective maturity ut- For sufficiently low values of ax 
the distribution is unimodal and centered around W = 1/2. Hence, in this early-time regime both options 
equally contribute to the total asset. However, if we allow the options to "mature" longer, we observe the 
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Figure 2: Exact distribution P{W), ©, for ax = 0.5 (dashed curve), ar = 1-63 (dotted), and ax = 3.5 
(dash-dotted). Thin solid lines are the corresponding asymptotic results in (IT9l ) (see Appendix IbJ. 

same surprising anomaly, which we have already encountered in the previous section: when ax reaches a 
critical value a^ ^ 1.63, the maximum at W = 1/2 ceases to exist and the distribution becomes close to 
uniform for 0.2 < W < 0.8. Thus any value of W in this range is nearly equally probable. For ax exceeding 
aj^, two maxima emerge continuously and the distribution changes its shape becoming M-shaped bimodal, 
with most probable values close to and 1. Although the average value W is still equal to 1/2 in this 
regime, the probability density now has a minimum at W = 1/2. Therefore, for /i = the path-dependence 
of the Asian-style variables does not suppress the transition to the disproportionate behavior but only shifts 
it to later times; q^ is more than three times larger than the corresponding value for the European-style 
variables. 



/i > 0. 

In this case the distribution function in ([8]l converges (Monthus and Comtet 1994, Comtet et al 1998), as 
ay — )• oo, to a limiting form defined by the first term of a discrete branch, 
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(10) 



Hence, for ax = oo the variable r has a very broad distribution characterized by a "fat" algebraic tail. For 
< /U < 1 this distribution is normalized but does not have moments. 

We do not attempt here to present an exact solution for P(>V) for arbitrary /i > and arbitrary ax- Instead, 
our aim will be to get a conceptual understanding whether P(>V) undergoes, at a certain unknown value of 
ax, a transition from a unimodal to a bimodal form. This question can be immediately answered if we find 
that the corresponding limiting form of PiyV) exhibits such a transition. 

Plugging the limiting form in (ITOl ) into (O and performing integrations, we find that in the limit a^ — )• oo, 
the probability density P(yV) converges to a limiting form 



P{W) ^ Poo(W) 
i.e., it tends to a beta-distribution. 
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The distribution on the right-hand-side of (fTTI ) has a different shape (modaUty) for < fi < 1, /i = 1 or 
fi > 1. For < // < 1 the distribution Poo(VV) is a [/-shaped bimodal. For /i = 1 it is uniform. Finally, 
Poo(W') is unimodal, centered at W = 1/2 for fi > I. 

Therefore, for < ^ < 1 the distribution P{W) will change its shape from a unimodal to a bimodal at a 
certain value a^. For the marginal /i = 1 case, the distribution P(W) will tend to a uniform distribution as 
aj- — )" oo. Finally, for /i > 1, the distribution will always remain unimodal and centered around the most 
probable value W = 1/2. Note that a similar transition was observed for the distribution of the occupation 
time on the positive axis for Sinai model with a drift jj. (Majumdar and Comtet 2002, Sabhapandit et 
al. 2006). Interestingly, this transition also occurs for jj. = 1. 



/i < 0. 

In this case the distribution ^(r) is given by the continuous branch in ([8]). Note that the expression in ^ 
can not be used directly, since it does not allow one to perform an integration over r in (|7]l. Thus, we will 
first try to obtain a plausible approximation for ^(r) valid for large ax, which will allow us to perform the 
integration over r. Then, on base of this result we will see whether the transition to the disproportionate 
behavior indeed takes place or not. After some straightforward manipulations (see Appendix IB]). we find 

where C is a constant (see AppendixiBl) and U{a, b, z) is the confluent hypergeometric function (Abramowitz 
and Stegun 1972). Note that for sufficiently large ax the approximation in (fT2)) works fairly well for any 
value of r (see Fig. |5]in Appendix |B]). Moreover, it exhibits exact asymptotic behaviors in the limits r — )• 

and r — )• oo. 

Equation (fT2l) predicts that ^(r) has a log-normal tail as r — ;• oo. Given such a slow decay, one expects 
that the integral in Q is dominated by large values of u. One finds then that, for sufficiently large ut, and 
W not too close to or 1, 

while for 1^ — ;• or 1 one finds that P{W) obeys Q (with possible logarithmic corrections, see Appendix 
IbI (fT9l)). Therefore, for /x < the distribution ^(lA^) does not reach a limiting form when ax -^ oo. This 
means, in turn, that /or all fi < 0, the distribution P(>V) exhibits a generic transition from a bell-shaped to 
an A/-shaped form, as ax passes through some critical value a^. 

In fact, this is a rather counter-intuitive result. The random variables At are integrals of a geometric Brow- 
nian motion - an exponential of a symmetric Brownian motion plus a constant drift term —a'^fi/2. One 
may naturally expect that for sufficiently large fi a contribution due to symmetric Brownian motion will be 
insignificant. This is precisely the behavior we observed in the /x > case, for which the transition to the 
disproportionate behavior takes place only for /x G [0, 1[ and is absent for /i > 1. Surprisingly, this is not 
the case for /x < and "disorder" embodied in Brownian terms Bt turns out to be always relevant, despite 
the fact that the constant drift term clearly provides a dominant contribution (for sufficiently large /x) in the 
exponential. 

We have performed numerical simulations of the race between two uncorrelated Asian-style random vari- 
ables [see (fT6l [TSl l below] which confirm our conclusions on the transition to a disproportionate behavior 




Figure 3: Probability density P(>V) obtained from numerical simulations (symbols) of the race between 
two uncorrected Asian-style variables for for fi = —1.7 and for two different maturities ax = 0.25 (a 
bell-shaped curve) and ax = 1-5 (an M-shaped curve). The critical value a^ ~ 1.12. The dashed lines 
correspond to a fit according to Q with ax replaced with a fitting parameter dx- 



for any // < 0. In Fig. [3]we plot P(>V) for jj, = —1.7 for three different values of the maturity ax, which 
clearly shows a transition from a unimodal to an M-shaped form, as the maturity exceeds a critical value 
a^ ~ 1.12. Guided by the above analytic argument (IT3] ). we have fitted these distributions using the expres- 
sion for the European-style variable, P{W) in Q, in which we replaced ax by some effective maturity dx 
used as a fitting parameter. As one may notice, the quality of the fit in Fig. [3] is very good. We also found 
numerically that a^ = a^ {fi) is a slowly decreasing function of ^u for ^u < 0, which implies that the larger, 
(by absolute value) // is, the earlier the transition takes place. 



3 Correlated increments. 



So far we have concentrated on the case of independent variables, to which one may object claiming that 
such a transition is spurious and correlations (Laloux et al 1999) between the increments dBl and dBl 
will "stabilize" the behavior of the European- and/or the Asian-style variables. In the remainder, we proceed 
to show that the transition to the bimodal shape is robust and the only effect of correlations is to shift the 
transition to later times. 

We focus on the case of very strongly correlated increments, when the penalty for having different dB^ 
and dBl grows with ut = dBi — dBi as uf. In this case the distribution function of the increments 



reads: 
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Physically, it corresponds to the situation of two Brownian particles coupled by a Hookean spring with 
rigidity l/x- The parameter x sets the scale of correlations. When x ~^ (^n infinite rigidity of the spring), 
the function F in ((14)) converges to a delta-function which signifies that the increments are forced to be 
equal to each other, dB\ = dB\ . Conversely, when x is large, the function F in (fT4l ) will tolerate large 
deviations of dB\ from dB\ . For x = oo we recover the Umit of independent variables. 

For F in (IT4l) . the distribution function P(VF) of the random variable W can be calculated exactly (see 
Appendix ICl): 

1/2 

and appears to have essentially the same form as the result for two uncorrelated European-style variables, 
^. The only difference is that now ut is renormalized by x. so that the transition to the bimodal form 
occurs at a|, = (2 + x^)/2x^ > 1/2. Hence, in presence of correlations the transition to the bimodal form 
is postponed for later times. 

For correlated Asian options, analytic calculations are much more involved and here we again resort to nu- 
merical analysis. To this purpose, we first generate two random walks Xn and Xn which evolve according 
to the following rules (in discrete time): 

^(1) _ ^(2) _ o 

4^^ = 4^1i - f -^ + ^i , 4^) = 4'2i - f -^ + 4^^ - > 1 , (16) 

where r/i'^ = dBf\ i = 1,2 and cr^ = E ([dBf\A = (1 + x^)/(2 + x^)- One has, of course. 



^(-)K|^)-"..^e..f-^.^f^,), 



E {r]n '%,' ) = if n 7^ n'. On the other hand, it is easy to see from (IT4l ) that the random variables 
Un = f]n — 7]n and Vn = vin + T^n are independent Gaussian random variables with correlations 



2x^ 
X^ 

(1) o„^ .(2) 



E{<) = -^ , E{vi) = 2 . (17) 
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Therefore, in order to generate the random variables ??„ and 1% - which in the limit of the large number of 
steps converge to dBi and dBl distributed according to (fT4l) . - we generate two independent Gaussian 
random variables u„ and Vn satisfying dTTl ). from which we obtain the desired ?]„ = (t>„ — Un)l2 and 

{'2) 

Vn = {vn + Un)/'i- Finally, from these two random walks it is straightforward to obtain two correlated 
Asian style options tj's from Q: 

" N 

r, = exp ^x« , (18) 

.n=0 



where N cxT. One can then compute the distribution P{W). By varying the parameters x and fi we obtain 
the phase diagram which is depicted in Fig. IH This is achieved as follows: 

• For /i > 1, we have shown previously that even in the absence of correlations, (which corresponds 
to X — ^ co), the distribution P(>V) is always unimodal so that it will remain unimodal for any finite 
value of X- In this case we have checked numerically that the distribution P{yV) converges to a 
A^-independent, unimodal curve for N in dTSl ) sufficiently large. 

• For < ^ < 1 and x ~^ oo, we have shown that the shape of P(W) changes from a bell-shaped 
to a [/-shaped form as A^ oc ot passes through a critical value Nc{fJ.,x ~^ oo), which we have 
numerically verified. For finite but sufficiently large x^ X > Xc{p) we continue to observe such a 
transition as N crosses a critical value Nc{n,x)' where Nc{ij-,x) increases as x is decreased, and 
eventually Nc{fi,x) diverges when x -^ Xc{p)- For x < Xc(/^) one observes instead that P{yV) 
converges to a A^-independent unimodal distribution. In our simulation, we thus identify Xc{lj) as the 
largest value of x for which, for a given < // < 1 fixed, such a convergence is observed. We have 
checked numerically that the "critical region" is very narrow, in the sense that Nc{n, x) grows very 
rapidly when x -^ Xc, which thus allows for a reliable estimate of the critical Une depicted as a dotted 
line in Fig. |4]in the (/i, 1/x) plane. Note, however, that a precise estimate of this critical line would 
certainly require a more careful analysis. 

• For /i < we observe instead a transition from an Af -shaped to a bell-shaped form of the distribution 
but the effect of correlations is qualitatively similar to the case < /i < 1. The critical line is 
identified as described above for < /x < 1. 
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Figure 4: A sketch of the phase diagram in the (/x, 1/x) plane. The dotted line starting at /x = 1 for 1/x = 
is the critical line ndx)- For fi > fic, the distribution P{yV) remains unimodal for all maturities. For 
fi < fic we depict the shape of P(W) for ax < a^, ot — ct^ and ot > a^- For < fi < ^c, the 
distribution P{yV) evolves from a bell-shaped curve to a [/-shaped one, while for /x < 0, P{yV) evolves 
from a bell-shaped curve to an M-shaped one as ax crosses a^. 



This phase diagram shown in Fig. |4]can be summarized as follows. For finite x, we find that P{W) always 
undergoes a transition from a unimodal to a bimodal form as the maturity passes through a critical value a^ 



for all fi < ^c = /ic(x)- As we expected, /Xc is a decreasing function of the strength of correlations: this 
is shown by the dotted line in Fig. |4] Similarly to the uncorrelated case, one also finds a different behavior 
for < fi < fic and /i < 0: in the first case, P{W) changes from a bell-shaped to a [/-shaped form, while 
f or // < it changes from a bell-shaped to an M-shaped form. This is depicted in Fig. IH Finally, for 
fi < /ic> we observe that a^ is an increasing function of the strength of correlations, similarly to the case of 
the European-style variables. 

We conclude with several words concerning a common feature of the transition to the disproportionate 
behavior observed in our paper. As one may notice, such a transition always takes place in situations when 
the first moment of the distribution function of the underlying S or A diverges as T — )• oo, and does not take 
place if the first moment remains finite at T = oo. It is a bit intriguing to see how the behavior at T = oo 
defines the transition which takes place at a finite T. 



A Independent European-style variables 

Let E {exp {—XW)}, A > 0, denote the moment generating function of the random variable W. The curly 
brackets here and henceforth denote averaging with respect to the distributions of Brownian motions B^ 

(2) 

and Bl . Explicitly, this function can be written down as 






Integrating over dBi, we formally change the integration variable Bi ^ W, which yields 

E \ e-^^ \ = , / — exp -XW In^ 

I / VMa^ Jo W{l-W) ^V 8ar VI -^ 

from which one immediately deduces Q. 



B Independent Asian-style variables 

The moment generating function of W can be written as 

E {e^^^} = / dTi dT2 ^(n) ^(ra) exp ( -A- 



n 



n +T2 

Integrating over dr2, we formally change the integration variable from T2 to W to give 
E L-^^l = f dW e^^^ / ri dn ^{Wn) ^ ((1 - W)ri) 
from which one may read off the result in ([7]). 
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/i > 

We note that in the limit W — )■ 1, P{yV) in ([6]) has the following asymptotic representation (Oshanin and 
Redner 2009): 

omn (l + ^^l^^) ( In^z tt^VWIuz , 

P(W) ^ V — ^ ^exp \ ^- ^ , (19) 



where z = {1 + \/>V)/-y/(1 — W). This asymptotic form agrees quite well with the exact result in Q, not 
only when W — )• 1, but also for moderate values of W (Fig.|2l). The asymptotics of V{yV) for W — )• can 
be obtained by merely changing W to 1 — lA^. 

It may be also worthy to remark that in the limit W ^ I the asymptotic form in ( |T9l ) follows, apart of a 
logarithmic factor In^'^ (1/(1 — VV)), the asymptotic form of the parental distribution ^(r) in ([D (Oshanin 
et al 1993b, Monthus and Comtet 1994, Comtet et al 1998): 

^(t) ~ — exp ( — In (r) I , r — )• oo 



For fixed W, the probabihty density P{yV) ~ l/y/or when ax -^ oo, which signifies that the large-T 
behavior of P{W) in (|T9l ) is supported by negative moments of r, which decay as £'(l/r") ~ l/y^or 
regardless of the order n (Oshanin et al 1993b, Monthus and Comtet 1994, Comtet et al 1998). 

/i < 

Whittaker function has the following integral representation: 



n-f + y" 
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= 2-,(-«./^e.p(-|)f^e.p(-^)K.„M, 

where Kiu{x) is the modified Bessel function (Abramowitz and Stegun 1972). Using the latter equation 
and a standard integral representation of Kj„(x), 



/•OD 

Kju(x) = / dt cos(nt) exp (— xcosh(t)) , 
Jo 



we find 



nr) = ^-'''"P(-°^^V4) 25^ ^,j^^_ ^,j^ ^^„, 



where r' = o^t j2 and il^iotT^ r) is given by 



,, , ^°^ dx ( a'^Tx^\.( aT\ ,^,, 
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with 

X, — ) = 3- / d^exp ( — xcosh {^) ) X (22) 

2 ^ (7raT)2 Jo V ar/ 

X vr cos — ^ sm 

V V «T / V «T / 

We now focus on ■(/'(aT,''") seeking for a plausible approximation in the limit a^ ^> 1. Note first that 
ip{aT, t) is a monotonically decreasing function of r and 

V'(aT,r = 0) = ^^^:^ + o(a^^/^). (23) 

In the large r limit, the integral in (|2TI ) is dominated by the small x behavior of the integrand. To obtain the 
small X behavior of 9 {x, ^) in the large ax limit, we take x = 1/(2 sinh {y^Jor)) and also change the 
integration variable ^ = u^fa^ in (l22l ). Setting then ay — >• 00 while keeping y fixed, we have 

X, ^) 2, -j^ye-y" = (24) 



3/2 



V2xy V "T \2xJ 



which defines the small x asymptotic behavior of 6{x, aT/2). 

On the other hand, for large x 

e{x,^)^2Ko{x)/V^a^'\ 

Hence, we may approximate [x, ^), for sufficiently large a-r, as 

Note that this approximate form reproduces correctly the exact behavior oiO {x,^) both for x — > (l24l ) 

and X — )• 00. 

We use next the small x asymptotic behavior of 9 (x, ^), (l24l ). to obtain from (1211 the following large-r 
asymptotic of ip{aT, t): 

The small-r asymptotic behavior of xp{aT,T) can be deduced from the large-x asymptotic of 9 (x, ^), 
which yields 

„„„.,.^!i^.W.t, (-1,1.4). (27, 

where U{a, b, z) is the confluent hypergeometric function (Abramowitz and Stegun 1972). 



12 



Combining the estimates in (l26l ) and (l27l) . we obtain the approximate form for ^(r) in (IT2l) with 



C 






■ exp 



ut^j- 



In Fig. [5j we compare our approximate expression for ^(r) in (IT2l ) against the numerical evaluations of 
the integrals in (l20l ) and (1211 for /i = —1 and different values of ar = 20, 30, 50. One notices that our 
approximation is quite accurate for all r > 0. As expected, the approximation works better as ax increases. 







200 



Figure 5: Probability density ^(r) for /i = — 1 and ax = 20,30,50 (from top to bottom). Each pair of 
curves have been rescaled by an appropriate factor so that they can be shown on a same graph. The symbols 
correspond to the exact value of ^(r), obtained by a numerical evaluation of the integrals in (l20l ) and (|2TI) . 
while the sohd lines correspond to the approximation given in (IT2l ). 



C Correlated increments 



We first divide the interval [0, T] into A^ subintervals dt approximating the Brownian trajectories B>p and 
B)p ' by the trajectories of random walks in discrete time (with time-step df) k: 



B. 



(n) 



N 



Y,dBf\n= 1,2, 



k=l 



jW 



where dB^ are the values of the increments at time moment k. In what follows, we will be interested by 
the behavior in the limit A^ — )■ cxd, dt — )• with Ndt = T kept fixed. 

Then, a random variable W, which defines a contribution of a given European-style option into the sum of 
two such options, can be written formally as 

M/ = ( 1 + n exp UdB'^^^ - dBJ^^)) ] . 



k=l 



13 



Then, for the generating function of this random variable we have 

N 



H^-}-(m) /:-/:n-r'-i^ 



(1)^r(2)x 



X exp 



<')' «')' «'-<'^" 



X exp 



2dt 2dt 2x^dt 



JV ._/^/^r(2)_^r(2) 



1 + nti exp [a [dB):> - dB): 

s)7:-/:n-i"-^ 



2 



xexp,-A!!fil-^,,-(i±^!M,x 






l + nfc=iexp(o-Ufc)/ 
which reduces, upon the integration over dB]^' , to 



^{-"'^^i;^£ /In-- 



X exp — exp 



2ff2di; ^\^ l + nf=iexp(cTnfc);' 
where g'^ = 2xV(2 + X^)- 

Further on, performing the integration over du^ we formally change the integration variable u^ — )• W, 
which yields 

1 JTT/ roo roo ^^1 



^ {-""} - ^rd:^ , mf^-^"- 1- L n "-- 



1 r dw 

a{27rg^dtf/^Jo W{l-Wy __ ■' ~^ k=i 

The latter equation implies that the distribution function P{W) in case of two European-style variables with 
correlated increments is given by 



14 



^^"p|-2i^(^^(V^)-E^'^^ 



Finally, by taking the advantage of the integral identity 



exp i-u^/2A) 1 f°° . 2 , N 

, = — / dbj exp (iuju — Auj 2) 

V2^ 2^J_^ ^^ ' > 



we readily perform the A^ — 1-fold integral in the latter equation to get the result in (fTSl) . 
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